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Quantum walk on distinguishable non-interacting many-particles and 

indistinguishable two-particle 

C. M. Chandrashekari'B and Th.Busch^'Q 
'Physics Department, University College Cork, Cork, Ireland 

We present an investigation of many-particle quantum walks in systems of non-interacting dis- 
tinguishable particles. Along with a redistribution of the many-particle density profile we show 
that the collective evolution of the many-particle system resembles the single-particle quantum walk 
evolution when the number of steps is greater than the number of particles in the system. For 
non-uniform initial states we show that the quantum walks can be effectively used to separate the 
basis states of the particle in position space and grouping like state together. We also discuss a 
two-particle quantum walk on a two-dimensional lattice and demonstrate an evolution leading to 
the localization of both particles at the center of the lattice. Finally we discuss the outcome of a 
quantum walk of two indistinguishable particles interacting at some point during the evolution. 



PACS numbers: 03.67.Lx, 03.67.Ac, 05.30.Ch 



I. INTRODUCTION 



The idea of a quantum walk, the quantum analog of 
the classical random walk, dates back to 1958 [l| and 
1965 [21 but the concept was formally developed only 
in 1990's 0-[al- In a one-dimensional situation a quan- 
tum walk evolving in position space spreads quadratically 
faster than its classical counterpart, due to the interfer- 
ence of amplitudes of the multiple paths 0, 0] ■ This was 
found to have interesting applications in quantum infor- 
mation theory, allowing for efficient quantum algorithms 
[8l4ll|. However, quantum walks have also been shown 
the be useful for coherent quantum control over atoms 
and quantum phase transitions [l^, to explain break- 
down phenomena in electric-field driven systems [13|, 
to give direct experimental evidence for wavelike energy 
transfer within photosynthctic systems |14l . [l5| , to gener- 
ate entanglement between two spatially separated system 
[Iq and to generate topological phases [I3|- Experimen- 
tal implementation of quantum walks have been reported 
using nuclear magnetic resonance (NMR) [iM[3i contin- 
uous tunneling of light fields through waveguide lattices 
[201 , the phase space of trapped ions [2l|, ^^ , single opti- 
cally trapped neutral atoms [2j] and single [2J, [23I and 
two-photon systems [20, [23 . All of these advances have 
made the area of quantum walks a very promising tool, 
just like its classical counterpart. 

Quantum walks are widely categorized into two forms, 
namely, continuous-time and discrete-time walks. In this 
article we will discuss the discrete-time quantum walk 
and in particular wc will focus on the quantum walk in a 
system of distinguishable particles. In Section [TTj we will 
define the distinguishable non-interacting many-particle 
quantum walk and discuss the dynamics and some of the 
results of the collective evolution of the many-particle 



system. We will show that for an evolution in which 
the number of steps is greater than the number of parti- 
cles, the collective probability distribution resembles the 
single-particle probability distribution. This can be ef- 
ficiently used for separating the different basis states of 
the many particle system and grouping them together in 
position space even when the initial states of the par- 
ticles are a randomized superposition of state. We also 
discuss the physical relevance of the study. In Section 
mil we look at the dynamics of a two-particle quantum 
walk and show that it is possible to localize the joint 
probability at the center of the lattice. We also look into 
the joint probability of the indistinguishable, both boson 
and fermion two-particle quantum walk evolution when 
the particle meet in the lattice after the walk evolution. 
We conclude in Section ITVl 

Though the review articles in this special issue intro- 
duce the concepts of quantum walks in great detail, we 
will briefly discuss the main features of the quantum walk 
evolution which will be relevant for this article here for 
self-consistency reasons. 

The discrete-time quantum walk of a single two-state 
particle in one-dimension is defined on a Hilbcrt space 
H = He ® "Hp, where He is the coin Hilbert space with 
the basis state described in terms of the internal state of 



the particle, | |) 
Hilbert space, Hn 



and I t> = 



The position 



Lp, has the basis states \ipj), where j G I 
is a set of integers associated with each lattice site. Each 
step in the evolution of the walk is described using a 
quantum coin operation 



B{9)^ 



_ ( cos(e') sin(6') 



sin(6l) -cos(6l) 



(1) 



which evolves the particle into a superposition of the in- 
ternal basis states and which is followed by the unitary 
shift operator 



^ ^ ^ [| ;)(; I ® |v^,-i)(v^,| + 1 t)(t I ® IV'.+iXv^j 
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(2) 



which transforms the state of the particle into a super- 
position in position space. Therefore, the full operation 
for each step of the quantum walk on the Hilbert space 
-He 



) T-Lp can be written in the form 



w{e) = s[B{e)®i]. 



(3) 



and the state after t steps of evolution is given by 

|*t)=W^(0)*|*i„). (4) 

Here 

|*i„) = (cos(V2)| ;) +e"'sin(<5/2)| t» ® IV'o), (5) 

is the initial state of the particle at a position j — Q. The 
coin parameter controls the variance of the probability 
distribution of the walk [y, 0, l^a [231 ^^.d the probability 
to find the particle at site j after t steps is given by 

P(j,t) = (^,|tr,(|vl/,)(*,|)|^,). 



II. DISTINGUISHABLE MANY-PARTICLE 
QUANTUM WALK 

Many particle quantum walks are fundamentally dif- 
ferent for systems of non-interacting distinguishable and 
indistinguishable particles. For the first one, the evolu- 
tion of the walk can be straight forwardly predicte d by 
considering many single- particle quantum walks [12, IS^l , 
whereas for the latter one many particle interference ef- 
fects, based on the bosonic or fermionic nature of the 
particles, strongly influence the evolution and makes it 
computationally hard to study [30, [31[ . 

Though the evolution of distinguishable particles does 
not involve many-particle interference effect, the collec- 
tive behavior of the single particle interference effects can 
reveal interesting features of the systems dynamics. Such 
systems can be approximately realized in cold, but ther- 
mal samples of neutral atomic gases in optical lattices 
[23| , which can be engineered to minimize the atom-atom 
interaction and dynamically control the atom transport 
|33l - l35| . Therefore, they have been sugg ested for observa- 
tion of quantum phase transitions |12| or for generation 
and control over spatial entanglement between different 
lattice sites |32l |. Another interesting question is the ex- 
ploration of the meeting probabilities and meting times 
of many-particles at pre-defined positions (see Ref. [3g 
for two particle meeting probabilities). 

In this section we will define the non-interacting distin- 
guishable many-particle quantum walk and discuss some 
of the interesting outcomes from the collective evolu- 
tion. To define a simple form of distinguishable many- 
particle quantum walk in one-dimension, we will con- 
sider a system of M non-interacting particles, where ini- 
tially exactly one particle occupies a lattice site and ev- 
ery particle has its own coin and position Hilbert space, 
n = {'Hc®Up)'^^ . If the number of particles is odd [3, 
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FIG. 1: (color online) Probability distribution of 51 particles 
initially located individually at positions j = —25 to j = -1-25 
after the quantum walk evolution using using the Hadamard 
operator Blji/'i) as the quantum coin, (a) The initial state of 
all the particles is -^=(1 \,)-\-i\ t)) and evolves in position space. 
The spread of the distribution with an increasing number of 
steps is clearly visible, (b) Initially all particles are either 
in state | 4-) or state | t) and subsequently subjected to the 
quantum walk of 500 steps. For all the particle in the initial 
state I 4-), the distribution with peak on the left is obtained 
and for all the particle in the initial state | X) , the distribution 
with peak on the right is obtained. 



the initial state can be written as 

■ _ A/-1 

;> + i|t) 



hP 



M 



V2 



•S) \-tJjj 



(6) 



which, after t steps, will evolve into 

;) + *it) 



*f> = [W{e) 



(^Mit 



V2 



•S) IV'j 



(7) 
Here W{9)^^ is the evolution operator for each step of 
the walk, which will evolve each particle into the superpo- 
sition of its neighboring positions, establishing the quan- 
tum correlation between the particle and the position 
space. After t steps, these correlations overlap resulting 
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FIG. 2: (color online) (a) The initial state has an antiferro- 
niagnetic ordering with neighboring particles being in differ- 
ent internal states, (b) Probability distribution of 51 parti- 
cles, initially in the state shown in (a) with one particle in 
each position ranging from j — —25 to j = -1-25, after a quan- 
tum walk of different number of steps using the Hadamard 
operator B{Tr/4) as the quantum coin. Lower parts in the left 
(right) of the distribution are due to the low contribution of 
state I i) ( I t)) from the particle initially in state | f) (I ^))- 
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FIG. 3: (color online) Probability of finding a particle in state 
j J,) of a single-particle walk with increasing number of steps 
if the particle was initially in state | \.}. For small values of 
the probability of finding the particle in state | \) is larger 
than for larger value of 6. 



[W{0) 



®M]ti^M 



[Aj-t\i)<E>\ijj-t) 



+Aj-t+i\ i) ® IV'j-t-n) + + A-+tl i) ® \^j+t) 

+B,-t\ t) ® \^j-t) + + B,+t)\ t) ® \^,+t)l (8) 



which can be written as 



j+t 



[W^W^"1*l*f^) « 



J2 [Ai\i)+Bi\t)]®\^. 

M-l \x=j-t 



(9) 



Here Al and Bl are the probability amplitudes of the 
state I I) and | t) of each of the particles initially at 
position j at the new position x, which range from (j — t) 
to {j + 1) after the t step walk. 

The probability distribution after t steps is given by 
the sum of the probabilities at a given lattice site of each 
particle k 



M 



Pij,t) = j2PkU,t) 



(10) 



k=l 



and the effective probability distribution for different 



numbers of steps is shown in Fig. 1(a) As expected, after 



t steps of the quantum walk, the M particles are spread 



between (M— i) and (M+t). In Fig. 1(b) the asymmetric 



probability distributions resulting from all particles be- 



ing initially in state 
From Fig. l(a)[ and 



l) or I t) after 500 steps are shown. 



1(b) it is clearly evident that when 



t ^ M the probability distribution profile of M particles 
resembles the single-particle profile. From earlier studies 
of single-particle quantum walks of t steps on a particle 
initially at position j = using Bg as the quantum coin 
it is known that the probability distribution spreads over 
the interval (— f cos(6'), tcos(6')) in position space and de- 
cays quickly outside this region [3, [2^ . For an M particle 
system the peak of the effective probability distribution 
is given by contributions from the probability of all M 
particles and therefore located at =F[tcos(0) — M/2]. 

Apart from all particles being initially in the symmet- 
ric superposition state of | \.) and | f) (see Eq.dS])), one 
can also consider a situation of antiferromagnetic order- 
ing (see Fig. 2(a)[ ), where two neighboring particles are in 
opposite states. In Fig. 2(b) we show the final probabil- 



ity distribution for this situation after a different number 
of steps. Though each particle undergoes an asymmetric 
evolution with the states | \.) moving left and the states 
I t) rnoving right, the collective distribution is symmetric 
due to equal number of particles initially in both states. 

Many-particle quantum walks of particles initially in 
antiferromagnetic order or in a completely randomized 
initial state are very useful for separating different ba- 
sis states of the particles in position space and grouping 
them together. Using a different angle 9 in the quantum 
coin operation one can find different outcomes for the 
probabilities of basis states grouped after the evolution 
of the quantum walk. To demonstrate this we will con- 
sider the examples of a single-particle initially in one of 
the basis state. 

A quantum walk of a single particle initially in state 
I \.) (or equivalently | t)) using a Hadamard coin {6 = 
7r/4) results in constructive interference towards the left 
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state after t steps can be written as 

t 

j=-t 

(11) 

where Aj.t and Bj,t arc given by the coupled iterative 
relations 



Aj,t = cos{0)Aj+i.t-i 
%t = -cos(0)S,_i,t_i4 



I- sin(6')6-,+i,t_i 
sin(6')^j_i,t_i. 



(12a) 
(12b) 



Straightforward algebra allows to decouple these equa- 
tions at the price of a time-dependence on the previous 
two steps 

Aj.t ^ cos{9) {Aj+i,t-i - Aj-i,t-i) - A,t-2 (13a) 
6,, t ^ cos(0) (6j+i,t_i - 6j-i,t_i) - Bj,t~2- (13b) 

By repeating this process of substitution one can find 
an expression linking Aj^t and Bj^t to the amplitude of 
the initial state of the particle and the angle, 9, of the 
coin operation. Therefore, the expression for the total 
probability of finding the particle in state | |) and | t) 
after time t is 
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FIG. 4: (color online) Probability distribution of 51 particles 
after 200 steps for different values of 6. The states | \) walk 
to the left and the states | t) to the right. In (a) the initial 
state of all particle was | J,) and in (b) the initial state of all 
particle was | t)- In (c) the initial state of each particle is 
randomly choose from | \,) and | '\). 



(right) of the origin and therefore localizes all particles 
with high probability on the left (right) of the origin. To 
understand this, let us look at the analytic form of the 
evolution after t steps using B{6) as coin operator. The 



To obtain a spatially symmetric probability distribu- 
tion for a particle initially in symmetric superposition 
state, the walk should be invariant under an exchange 
of |0) -O- |1), and hence should evolve Aj^t and Bj^t alike 
(as, for example, the Hadamard walk does [33 )■ From 
the above analysis we see that Aj^t and Bj^t are sym- 
metric to each other and evolve alike for all value of 9 
only when the initial state of the particle is a symmetric 
superposition state. When the initial state is | \), the 
walk will evolves with constructive interference towards 
left and destructive interference to the right, (exact form 
depending on the value of 9) and vice versa when the ini- 
tial state is \^). The associated probability amplitudes 
oscillate strongly between the left and the right hand side 
for small numbers of steps and stabilize for longer times. 
This can be seen in Fig. [3] and also directly from Eq. p^ 
when realizing that the amplitude at each position oscil- 
lates and the range of oscillation reduces as the amplitude 
at each position decreases over time [43 |. 

For a particle initially in state | \) a smaller value of 9 
returns a high probability of finding the particle in state 
I I) but if the initial state is | f) the probability of finding 
the particle in | |) will be very low. We should also note 
that a small probability of state | t) (I D) is present 
along with the state | |) (| f)) to the left (right) of the 
origin but that will not alter the trend. From the above 
analysis we can conclude that an initially randomized 
many-particle state can be efficiently sorted in position 



space with respect to its basis states. This in turn aUows 
to create an ordered state with high probabihty. 



To demonstrate this we show in Fig. 4(a) the proba- 



bihty distribution for different values of 6 for a sample 
of 51 particles after 200 steps when the initial state of 
all the particles was | \) and Fig. |4(b)] shows the same 
for an initial state of | f)- A strong asymmetry is visi- 
ble for both cases. In contrast, the probability distribu- 
tion shown in Fig. 4(c) assumes that the initial state of 



each particle was randomly chosen from | J,) and | f) and 
the anisotropy in the final distribution vanishes. From 
Figs. 4(a) and ( |4(b)[ ) one can also see that for increasing 
9 the probability distribution widens and its maximum 
amplitude decreases. 
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FIG. 5: (color online) Probability distribution of two distin- 
guishable particles on a two-dimensional lattice using Bin /A) 
as quantum coin operation, (a) Joint probability distribution 
of two particles staring at (0, 0) and (20, 20) with initial states 
I \) after 10 steps, (b) Localization of two-particle probabil- 
ity distribution at the center of the lattice after a one time 
bit-flip operation on both particles at i = j/2 was introduced. 




(c) 



FIG. 6: (color online) Quantum walk on a two-dimensional 
lattice for two indistinguishable bosons initially at (0, 0) and 
(20, 20) in I 4-4-) a-nd interacting via ax after 20 steps using 
B{-n /A) as the quantum coin operation. The probabilities for 
finding the particles in the state (a) | 4-4) (b) | tt) ^'iid (c) 
I 4t) a^re shown at i = 20. The relative height of the final 
distributions however can be shown to depend on the initial 
state. 



III. JOINT PROBABILITY OF TWO-PARTICLE 
QUANTUM WALK 

Two-particle quantum walks have been studied from 
various perspectives [38l44l| and first experimental im- 
plementations have recently been reported |22l . |26| . Here 



we will discuss the probability distribution of a quan- 
tum walk using two distinguishable particles on a two- 
dimensional lattice and present a protocol to increase the 
meeting probability of the two particle at a particular lat- 
tice after a particular time. We then compare this to the 



quantum walk evolution of two indistinguishable parti- 
cles which only interact at the end of a certain number 
of steps. 



Each step of the evolution can be implemented by 
W2{e) = [B{e) ® 5*1] ® [B{9) ® 52] and after t steps the 
state is given by 




FIG. 7: (color online) Probability distribution of finding the 
two fermions staring at (0, 0) and (20, 20) with initial states 
I 4,) and interacting with a^ after 20 steps of quantum walk on 
a two-dimensional lattice using B{tt/4) as quantum coin op- 
eration. The state | 4-t) is the only possible state for fermions. 

To define a two-particle quantum walk we will consider 
a two-dimensional square lattice and label the two axis 
as X and Y such that [x, y) represent a position on the 
lattice. We will consider two particles initially in state 
I ^) at diagonally opposite points (0,0) and (j,j), 



I* 



= [\\)® IV-O.O)] ®[\i)® \4>0,3 



(15) 



where j is the length of the lattice, which has j x j posi- 
tions. The shift operator for the quantum walk evolution 
is defined separately for both particles in such a way that 
they evolve towards each other, 

x,y 
+ lt>(tl®l^x,y+l)(V'x,y|] 



^,y 



[w^2(^)]*|*L> = {[B{e) ® ^1] ® [B{d) ® 52]}* 

x{[U>®k'o,o)]®[|i)®|V'w)]}- (17) 



The two particles meet each other for the first time after 
t ~ j steps and the meeting probability at each position 
is different for the distinguishable and indistinguishable 
case. 

Two distinguishable particles: In this case the joint 
probability of the two particles at each position at the 
time of meeting each other is the sum of the probabilities 
of both individual particle. In Fig. 
distribution for both particles after t 
on a 20 X 20 lattice. The first time the two distributions 
overlap is at i = 20 where they spread along the diagonal 
of the lattice (not shown). If, however, we introduce a 
one time bit-fiip operation, ax at t=j/2=10, 



5(a) 



we show this 
7/2 = 10 steps 



[w{e) 



ig,2]]/2r 



[w{0y 



i2lj/2\^2 \ 

\ I 1713/ 



(18) 



+lt>(tl®l^.,J,-l>(^.,J 



(16) 



one can see from Fig. |5(b) that the evolution can be 
reversed, which leads to localization of both the particles 
in the center of the lattice at time t = j with a good 
probability. 

Two indistinguishable particle: If the two particles are 
indistinguishable, their probability distributions interfere 
when they overlap at the same position in the lattice. For 
bosons the allowed states at each position in the lattice 
are | H) , \ ft) or | It) = | |t) i whereas for fermions these 
are restricted to | |t) = | |t) • The probabilities for these 
states to be obtained at each position at time t are then 
given for bosons as 






Pi 



n) 



1-^^*1' -14*1' 



EMW ■ 14*1' + Ht^ ■ \^U' + \AU' ■ \Bl\'+Al\^ ■ \Bl,\^] 



Ht\'-\Bl\' 



itt) E,[IA-*I' • 14*1' + Ktl" ■ \BU' + 14*1' • l^'tl' + 4*1' • l^x*l'] 



14*1' -l^'*!' 



14*1' -l^i,*!' 



E,[I4*I' • 14*1' + Ht\' ■ 14*1' + 14*1' ■ 14*1' + 4*1' ■ 14*1'] ■ 

I 



(19a) 
(19b) 
(19c) 



Here A'^ and A'' are the amplitudes of the particles a and b to be in state | |) , and S" and B^ are the amplitudes 



to be in the state | f)- We show these probabiHties two 
particles initially at (0,0) and (20,20) and meeting after 
evolving for 20 steps of walk using Eq. (fT7|) in Fig. 6(a)[ (b) 
and (c). 

If the particles are fermions the probability of finding 
the two-particle in the only possible state at each posi- 
tions is 






\Al,\'-\BU-' + \A 



j,t\ 



\Kt\ 



E. 



\\Al,\^-\Bl,\^ 



\A) 



\? 



\Blt?] 



(20) 



which is shown in Fig. [7] for the same parameters as in 
the bosonic case above. The difference to the bosonic 
case is clearly visible. Using different initial states of the 
particle or different coin operations during the evolution 
will of course alter the probability distribution. Intro- 
ducing a one time bit-flip operation half way through the 
evolution for indistinguishable particles as we did for dis- 
tinguishable particle will lead to localization of the join 
probability at the center (not shown). 

From this one can see that even a one time particle- 
particle interaction in an indistinguishable many-particle 
quantum walk can result in different probability distribu- 
tions which might be useful for applications in quantum 
information and other fundamental quantum mechani- 
cal experiments. With the possibility of increasing the 
number of steps, the number of particles and the num- 
ber of time the particle-particle interaction is introduced, 
the evolution gets even more interesting and complicated, 
but becomes computationally difficult. Recently,for the 
case of two atoms in an optical lattice performing a quan- 



tum walk with interactions via cold collisions the appear- 
ance of a bound state has been predicted ^^ , which gives 
scope for further exploration of the dynamics using our 
approach for many-particle system by introducing inter- 
actions at regular intervals. 



IV. CONCLUSION 

We have presented a number examples of quantum 
walk dynamics of many-particle system with different 
initial states of the particles. Though the distinguish- 
able many-particle quantum walk dynamics does not in- 
volve many-particle interference during the evolution we 
have shown that it can be effectively used to separate 
the eigenstates of the particles position space and group 
them together. We have also presented an example of 
two-particle quantum walk dynamics with defined inter- 
action that can lead to localization of two distinguish- 
able particles at the center if they start their walk from 
opposite ends of the lattice. Extending this scheme to 
indistinguishable boson and fermion pairs results in the 
different probabilities for finding the two particles in the 
allowed combination of states. Recent experimental de- 
velopments in implementing quantum walks and using 
quantum walk models to simulate and understand some 
of the dynamics process in nature suggests that collec- 
tive dynamics of many-particle system will very useful 
for further studies. 
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